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Abstract. Microtearing modes (MTMs) are unstable in the shallow gradient 
region just inside the top of the pedestal in the spherical tokamak experiment 
MAST, and may play an important role in the pedestal evolution. The linear 
properties of these instabilities are compared with MTMs deeper inside the core, 
and further detailed investigations expose the basic drive mechanism. MTMs 
near the MAST pedestal top are not well described by existing theories. In 
particular the growth rate of the dominant edge MTM does not peak at a finite 
collision frequency, as frequently reported for MTMs further into the core. Our 
study suggests that the edge MTM is driven by a collisionless trapped particle 
mechanism that is sensitive to magnetic drifts. This drive is enhanced in the outer 
region of MAST by a high magnetic shear and a high trapped particle fraction. 
Observations of similar modes on conventional aspect ratio devices suggests this 
drive mechanism may be somewhat ubiquitous towards the edge of current day 
and future hot tokamaks. 



1. Introduction 

Initial analytic studies suggested that tearing modes should be stable at high binomial 
perpendicular wavenumber, k y , due to increased field line bending [1], leading to a 
focus on larger scale, "gross" tearing modes. A kinetic study of the tearing mode 
found that an energy dependent collision operator could lead to an additional drive 
from the electron temperature gradient [2]. This drive can overcome the stabilising 
influence at large k y , allowing unstable microtearing modes (MTMs) to exist. The 
parameter v = v e i/oj, is important for MTMs, where v e i is the electron-ion collision 
frequency and uj is the frequency associated with the mode. Analytic treatments in 
the collisional (v 1), semi-collisional (p > 1 and v <C k^v^ h e /u> 2 ) and collisionless 
(v « 1) regimes were performed in slab geometry [3]. Here the instability arises due to 
the time dependent parallel thermal force providing an asymmetry in the parallel force 
on electrons as a consequence of the energy dependence of the collision operator [4, 5] . 
In the absence of collisions this asymmetry disappears and the slab drive is therefore 
expected to vanish for sufficiently small v. Likewise large v becomes stabilising as the 
collisions prevent the electrons from building a perturbed current. Kinetic calculations 
in large aspect ratio toroidal geometry revealed an additional drive mechanism that 
depends on trapped particles and is effective in the banana regime, v < r/R = e [6]. 
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The trapped particles themselves do not carry the perturbed current but collisions 
generate the instability by allowing current to grow in the barely passing particles close 
to the trappcd-passing boundary, in a process that remains effective even for 
This effect is found, for realistic tokamak parameters, to be lost at higher collision 
frequencies in the collisionless regime, v < 1 < v/e [7]. Importantly combining the 
slab and trapped particle drive mechanisms leads to an MTM growth rate which 
peaks for v ~ O (1). Observations of a strong inverse collisionality dependence of the 
thermal confinement time made on both MAST [8] and NSTX [!)] may be consistent 
with simulations showing MTM driven transport increasing with v e i [10]. 

Fully electromagnetic gyrokinetic simulations are now able to study microtearing 
modes numerically in experimentally relevant scenarios. Linear gyrokinetic 
calculations have found unstable MTMs in a wide range of equilibria including at 
mid-radius in spherical tokamaks (STs) [11—16], in simple large aspect ratio shifted 
circle model equilibria [15] and towards the edge in ASDEX Upgrade [17]. MTMs 
exhibit tearing parity, where in ballooning space the perturbed parallel magnetic vector 
potential, A|| , is even in the ballooning angle 9. This is associated with reconnection 
at the rational surfaces generating small scale island structures. 

When the amplitude of these islands is sufficient they will overlap to generate a 
stochastic field, which gives rise to significant electron heat transport [18]. Estimates of 
the electron thermal diffusivity in NSTX based on a model of stochastic field transport 
[19] are found to be within a factor 2 of the experimental levels over a region in which 
MTMs are the dominant instability [20]. The first successful nonlinear simulations 
of microtearing turbulence [21, 22] indicate that, in the absence of sheared flows, 
the associated electron heat flux can indeed be significant, and within the range of 
experimental observations. The effect of sheared equilibrium flows is not yet clear, 
with conflicting findings emerging from these two studies. The dependence of this 
microtearing turbulence on electron beta, j3 e — 2fion e T e / B 2 , inverse temperature 
length scale, l/£r e = T' e /T e , and collision frequency, v ei [10, 23] is broadly consistent 
with previous linear studies [15—17] . Whilst these gyrokinetic simulations agree 
qualitatively with the two drive mechanisms discussed earlier, through the dependence 
of the MTM growth rate, 7 MTM , on i>, the existence of a critical T' e and the observation 
that u ~ w» e , there is evidence that magnetic drifts, which have not been adequately 
treated analytically, are also important [15]. In particular, the energy dependence of 
the collision operator is vital for both analytic drives but in numerical simulations 
this had little impact in the presence of magnetic drifts [15]. Indeed it was found that 
both magnetic drifts and perturbed electrostatic potential, 4>, could be destabilising, 
and in the absence of both of these effects the MTM was found to be stable [15]. 
It is likely that there are multiple mechanisms occurring simultaneously to drive (or 
damp) MTMs, with the local parameters determining the relative contribution of 
each mechanism. This can lead to different scaling of 7 MTM , depending upon which 
mechanism is dominant. For example, Ref. [23] notes that <p is destabilising for 
low safety factor, q < 3, but stabilising for q > 3, suggesting that the dominant 
driving mechanism may be undergoing a transition as q varies. Whilst two MTM 
drive mechanisms have been uncovered by analytic theory, it is clear that additional 
mechanisms, related to the magnetic drifts, are absent from the existing literature. 

Recent linear gyrokinetic studies of the pedestal in MAST [24] and JET [25] 
utilising the fully electromagnetic initial value gyrokinetic code GS2 [26] have found 
unstable MTMs in the shallow gradient region at the top of the pedestal, where they 
may play an important role in the pedestal evolution [27, 28]. These edge MTMs 
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exhibit some significant differences to the more famiiiar core MTMs, and this paper 
provides an in-depth study of the edge modes. A clear difference can be seen in (/>, 
which is much less extended in 9 in the edge case than in the core, as shown in figure 1. 
In both cases A|| is strongly peaked about 6 = 0, decaying by 6 = ±7r. This is a more 
extreme example of the comparison of MTM cigenfunctions on surfaces at r/a = 0.6 
and r/a = 0.8 reported in Ref. [16]. It should be noted that the magnetic shear, s, is 
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Figure 1. The electrostatic potential perturbation, 0, due to MTMs from linear 
GS2 simulations of MAST at (a) i/>jv=0.58 (shot #27905) and (b) ip N =0.94 
(composite based on shots #24452, #24459 and #24763, see Ref. [21]). 



much larger in the edge, and both cigenfunctions extend in the ballooning coordinate, 
6, to approximately the same normalised radial wavenumber, k x pi = §6k y pi ~ 140, 
where pi is the ion Larmor radius. This is consistent with the mode extending until 
k x 5o ~ O (1), where Sq is the semi-collisional width [3]: 

s = Ls y^*^ (i) 

with the shear length, L s = Rq/s. These cases are both in a similar collisionality 
regime as v = 0.26 and 0.44 for the core and edge cases respectively. 

In section 2 we introduce a local equilibrium from the MAST edge that is unstable 
to MTMs, and reduce this to a simpler model equilibrium with similar microstability 
properties. This provides a reference equilibrium for detailed linear gyrokinetic studies, 
presented in section 3, that probe the basic driving mechanisms for MTMs at the edge 
of MAST. Final conclusions are presented in section 4 



2. Equilibrium parameters and simplifications 

We base our studies on a reference local equilibrium from the plateau region at the top 
of a MAST H-mode pedestal, which is unstable to MTMs \. The reference flux surface 
is tpN = 0.94 at the midpoint during the ELM cycle, with the equilibrium parameters 
given in table 1. The growth rate spectrum is shown in figure 2(a). The minimal 
equilibrium conditions necessary to drive MTMs unstable are sought by progressively 
simplifying the equilibrium assumptions. 

Sensitivity to flux surface shaping is investigated by fitting the reference MAST 
equilibrium using the a- a shifted concentric circle model [29]. This model allows easy 
independent control of the geometrical parameters affecting trapping, e = r/R, and 

\ A full account of the equilibrium reconstruction from MAST data is given in Ref. [24]. 
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Table 1. Equilibrium parameters characteristic of MAST shot 24763 at the mid- 
point in time between two ELMs for i/ijy = 0.94. 

q s r/R L ref /L T L r /L n e a v ei (v th4 /L ref ) 

4.66 7.67 0.805 5.88 0.36 0.015 -5.66 1.98 



the magnetic drifts, R, which shall later be shown to be crucial to the MTM instability 
studied here. The circle is a crude fit to the edge of MAST, as illustrated in figure 2(c) 
which compares this fit with the experimental flux surface. The 7 M tm spectrum for 
the circular fit is shown in figure 2(b). There is a significant shift in the k y pi at which 
7mtm peaks relative to the shaped surface case, but the magnitudes of 7 MTM (and ui) 
are within a factor 2. This is consistent with previous studies showing that shaping is 
not essential for MTMs [15]. 




en 



do 




Figure 2. (a) 7mtm spectrum for the MAST flux surface at ipw = 0.94, 
highlighting the peak wavenumber (□) corresponding to the <f> eigenfunction of 
figure 1(b). (b) 7mtm spectra for the circular equilibrium fit [- — in (c)] with: 

the standard full physics model excluding By ( ); adiabatic ions (*); including 

-By (□) and neglecting <f> ( ). (c) Circular fit ( ) to the tpN = 0.94 flux surface 

( ) along with the last closed flux surface ( — ■ — ). 



Figure 2(b) shows that calculations with fully kinetic and purely Boltzmann ion 
responses yield very similar 7 M tm spectra. Previous simulations of MTMs in the core 
also found that 7mtm is insensitive to including fully kinetic ions as the ion response 
is close to Boltzmann [15-17]. In early kinetic treatment of the collisionless and semi- 
collisional MTMs the ion response was neglected [.3], but can be strongly stabilising 
when pi> d [31 )] , where d is the width of the current layer associated with the mode 
(e.g. the semi-collisional width, 6o). Calculating the semi-collisional width for our 
reference equilibrium from (1) gives o~o/Pi ~ O (lO -3 ) and with this current layer 
width the stabilising impact of the ion response is expected to be important [30] . 

The growth rate is insensitive to including comprcssional magnetic perturbations, 
Bn, and only weakly sensitive to including the electrostatic potential, (j>, which is 
slightly stabilising and has a modest impact on the frequency spectrum (not shown). 
Subsequent simulations in this paper will use the s-a equilibrium model, neglect the 
nonadiabatic ion response and B\\, and retain <fi due to its impact on ui at large k y pi. 

The studies of section 3 are based on scans around the reference equilibium, during 
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which MTMs can become subdominant to other instabilities. GS2 is an initial value 
code, and subdominant MTMs are tracked, with this up-down symmetric equilibrium, 
by filtering to keep only the component of the nonadiabatic perturbed distribution 
function with odd parity in the parallel direction. These are tearing parity modes (i.e. 
modes where </> is odd and A|| is even about 9 = 0). Finally both the semi-collisional 
width, 5q, and the collisionless width, S n = p e y^2/ /3 e , are resolved by using a domain 
that is sufficiently extended in 9, (— lln < 9 < lln). 

3. Linear mode analysis 

There have been several linear gyrokinetic studies of how MTM stability depends on 
equilibrium parameters [15-17]. Here we explore a new region of parameter space, by 
moving to extremely low aspect ratio and high magnetic shear, which are typical of 
the MAST edge. 

3.1. Temperature and density dependence 

A finite electron temperature gradient is essential for both MTM drive mechanisms 
described in section 1, with the onset of instability arising above a certain threshold. 
The MTM real frequency, w, is predicted to vary approximately linearly with the 
electron diamagnctic frequency, w* e , with a precise relationship that depends on the 
driving mechanism. 

Scans have been performed by varying the normalised gradient length scales, 
L re f/LT e and L re f/L nc (with L re t some reference length scale and L„ e = n' e /n e ), 
independently, at fixed values of all other parameters §. The resulting growth rate 
spectrum in figure 3(a) show a finite threshold temperature gradient that increases 
with k y , and that 7 MTM increases monotonically above this threshold. Figure 3(b) 
shows that MTMs are unstable at L re f/L 7lc = 0, and that 7 M tm is maximised at 
finite L re f/L ne similar to previous findings [15, 16]. In both scans u is found to 
be reasonably well described by w* e (a + brj e ) where rj e = L n< ,/Lx a , in qualitative 
agreement with analytic predictions. 




(a) (b) 



Figure 3. The growth rate as a function of k y pi for varying (a) temperature 
and (b) density gradient length scales. The location of marginality, 7mtm = 0, is 
given ( ). 

§ The normalised pressure gradient, a = Rq 2 /3/L p , was held constant in these scans. 
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3.2. Beta 

7mtm is insensitive to 0, and MTMs are driven by the magnetic perturbation, 
A||. Therefore 7 M tm will be strongly affected by /?, which controls the strength 
of magnetic perturbations through Amperes law. Electromagnetic instabilities, like 
kinetic ballooning modes and MTMs, are typically unstable above a threshold j3, with 
growth rates that then increase strongly with (3 [31]. This may explain discrepancies 
between different tokamaks in the observed confinement scaling with /3 [32]: increases 
in j3 that cross the threshold will increase transport whilst increases that remain below 
the threshold will have less impact (and may stabilise other instabilities [33, 34]). 

The growth rates are shown for a range of k y values in figure 4(a) for varying 
/? e (scaling a consistently). The variation seen here is somewhat complex, though 
some general remarks can be made. There is a clear stability threshold in (3 e , which 
increases approximately linearly with k y , above which 7 MTM increases rapidly with (3 e . 
For sufficiently high /? e , further increases in (3 e become stabilising, as was also seen 
in Ref. [15]. Varying a consistently with j3 e significantly enhances the stabilisation 
at higher /3 e , which is consistent with magnetic drifts becoming more favourable at 
higher a [35] . Holding a fixed in this scan, 7 MTM decreases more slowly with increasing 
/3 e above the most unstable value, and furthermore the instability threshold moves to 
lower j3 e . The local minimum in 7 MT m at j3 e ~ 0.022 is only seen in the scan with a 
varying consistently, and not at fixed a. 

kyp i= i.Q 
kyp i= o.9 
k v p i= 0.S 
k y pi=0.7 
k v Pi=Q.Q 
k y p i= o.5 

kyp i= 0A 
kyp i= Q2 
k y pi=0.2 

0.01 0.02 0.03 lxlCrtxld 1 1x10° lxlO 1 lxlO 2 lxlO 3 

Pe V ei (vth,i/L re f) 
(a) (b) 




Figure 4. 7mtm as a function of k y pi and (a) /3 e and (b) v £ i. During the (3 e 
scan the normalised pressure gradient, a, is varied consistently. The results of 
simulations with u e i = are shown in (b) by the filled semi-circles on the y-axis, 
indicating a substantial growth rate even in the absence of collisions. The key to 
the kypi values also applies to figures 5, 6, 7. 



3. 3. Collision frequency 

It has already been pointed out that collisions play an essential role in the existing 
analytic drive mechanisms for MTMs. Linear gyrokinetic simulations have generally 
reported growth rates that peak for v ~ O (1), as may be expected from a mode driven 
by a combination of slab and trapped particle drives. Recent simulations have shown 
7mtm dropping by only a factor of 2, as v e i falls by over two orders of magnitude from 
its value at the peak [23] , which suggests that as the collision based drive is removed, 
a further substantial drive mechanism remains. 
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Fi gure 4(b) shows 7mtm &s a function of v e i for each ky. Increasing the collision 
frequency well above v ~ 0(1) is stabilising. It is more striking that 7 MT m for the 
dominant mode, and at several other values of k y , does not peak at finite is e i, but 
remains constant or even slowly increases as v e % decreases all the way to zero, which 
is in stark contrast to the "usual" behaviour ||. The trapped particle drive mechanism 
of Rcf. [(>] must vanish at v e i = as it requires collisions, and cannot therefore be 
responsible for the instability seen here. A collisionless mechanism is required, which 
cannot rely on the time-dependent thermal force. 

Ref. [15] reported 7 M tm peaking at v ~ O (1) for mid-radius local equilibria in 
MAST, but closer to the MAST edge we find that 7 mt m peaks at v ei ~ 0. Could 
this be explained by significant differences between the values of inverse aspect ratio, 
e = r/R, and magnetic shear, s, at these two locations? 

3.4- Aspect ratio (trapped particles) 

Varying only the inverse aspect ratio, e = r/R, in the s-a model, corresponds to 
changing the trapped particle fraction whilst holding all other parameters fixed. This 
effectively scans in minor radius, r, at fixed R. The results from this scan, illustrated 
in figure 5(a), reveal a strong dependence of 7mtm on e, and suggest that trapped 
particles are important to the linear drive for the reference value of v e i- The decline 
in 7mtm with decreasing e is nearly uniform for k y > 0.5, but rather weaker at lower 
k y . This suggests that at low k y the trapped particle drive may be complemented by 
another mechanism at the nominal v e i, which would also be consistent with figure 4(b). 

Ref. [15] found that trapped particles are destabilising for v <C 1, but stabilising 
for v > 0.5 which is the reference collisionality regime here. Furthermore, the trapped 
particle drive was shown to be most effective at low e (unlike in figure 5(a)). This is 
consistent with the trapped-passing boundary providing an instability drive, whilst the 
trapped particles are themselves stabilising owing to their inability to carry the current 
perturbation. The situation is different for MTMs at the MAST edge. Figure 5(a) 
suggests that trapped particles provide a direct MTM drive, and figure 4(b) shows 
that this survives without collisions. 




6 Vei (Vth,i/L re f) 

(a) (b) 



Figure 5. 7mtm as the inverse aspect ratio, e = r/R, is varied for (a) v e i = 1.98 
at each k y (k y key given in figure 4) and (b) k y pi = 0.6 over a range of i/ e i, with 
— indicating marginal stability. The frequency, ui, remains between -1.5 and 
-2.8 throughout these scans. 

|| 7mtm for the lowest k y mode peaks at finite v e i, and drops only slowly with decreasing v e i, closely 
resembling the dependence presented in Ref. [23]. 
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Fi gure 5(b) shows how ^ MTM for the dominant mode, at k y pi = 0.6, depends on 
v e i and e. The 7mtm dependence on e is strongest at low v e i. At low e the growth rate 
maximises at v e i ~ O (10), but at high e the growth rate peaks at the minimum v e i. 
Indeed for e < 0.3 the MTMs become stable for sufficiently small v e i, and a strong 
peak in 7 MTM is seen for v ~ 10, which is consistent with previous findings [15, 16, 23]. 
The role of e in enhancing 7 MTM should be most important near the pedestal region 
where the trapped particle fraction variation is maximised, especially in STs where e 
approaches 1. 

3. 5. Safety factor and magnetic shear 

The safety factor, q, and magnetic shear, s, are significantly larger in the edge plateau 
of MAST than at the mid-radius surface studied in Ref. [15]. Figure 6(a) illustrates, 
for a range of k y values, the complicated dependence of 7 MTM on q. 7 M tm exhibits 
multiple peaks at different q values, and these positions vary with k y : e.g. the k y pi=l 
mode has growth rate peaks at q ~ 3 and q ~ 9, and a local minimum at q ~ 6. 
Interestingly the MTM is stable for q < 2. These features may be related to the 
impact of q on bounce and transit frequencies, which are inversely proportional to 
q. Figure 6(b) shows results from a scan in s, and clearly indicates a preferred value 
of s for which 7 M tm is maximised for each mode. The most unstable s decreases as 
k y increases, and at lower s the peak of the 7mtm spectrum moves to higher k y pi. 
Changes in s affect the magnetic drift frequency, ujd, which in the next section will 
be shown to impact on the growth rate. 




2468 456789 
q I 

(a) (b) 



Figure 6. 7mtm as a function of (a) q and (b) s. The k y key is given in figure 4 



3. 6. Drift frequency 

In the s-a model the curvature and V£? drifts are equal, and for zero ballooning angle, 
6 = 0, are proportional to 

D (9) = ei [cos (0) - (a sin (0) - s0) sin (0)] (2) 

where e; = 2L re f / R determines the strength of the drifts. A scan in drift frequency 
was performed by varying e; around its reference value e; = 1.435 with all other 
parameters fixed. Figure 7(a) shows that the growth rate peaks at a particular q, 
which varies with k y , and that 7 MT m is more sensitive to the drifts at higher k y . 
The peak growth rate occurs at a drift strength factor that decreases approximately 
linearly with k y , suggesting an optimal value of u>d for which 7 MTM of each mode is 
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maximised. This indicates that some form of drift resonance may be important. In 
the absence of magnetic drifts (i.e. e; = 0) the MTMs are stable, showing that the 
slab drive is insufficient for instability 

Independent scans in the magnetic drift frequencies for trapped and passing 
particles, reveal that 7 MTM is only sensitive to the trapped particle drifts and that 
passing particle drifts arc unimportant. We note that the previous scan in s was 
effectively a scan in the radial component of u>d, as this is oc sm(9)§9 (from (2)). 
Therefore the similarity of figures 6(b) and 7(a) indicate that the radial component of 
the magnetic drift is the dominant influence on the drive mechanism. In ballooning 
space the radial wavenumber exceeds k y for §6 > 1, which arises for 9 > 0.13 in the 
edge, and for 9 > 3.49 for the mid-radius MAST parameters of Ref. [15]. The radial 
component of the drift frequency for trapped particles is clearly more significant at 
the edge of MAST than at mid-radius. Trapped particles and their radial drifts seem 
to play an essential role in the MTM drive mechanism at large e. Analytic theories 
of the MTM either neglect the magnetic drift frequency, u>u, ot assume ojd w. For 
MTMs at the edge of MAST it is demonstrated here that the magnetic drifts cannot 
be neglected. 




0.5 1.0 1.5 2.0 2.5 1 2 3 4 

e« V±/V t h,e 
(a) (b) 



Figure 7. (a) 7mtm as the normalised major radius, e; = 2L re f/R, is varied for 
fixed e (k v key given in figure 4) and (b) the location in velocity space where u 

for kyPi = 0.6 matches the drift (evaluated at 8 = 0, ), precession ( - ■ — ) 

and bounce ( ) frequencies. 



3. 7. Frequencies 

It is of direct interest to analytic theory to ask how the MTM mode frequency, u, 
compares with the natural trapped particle frequencies: the bounce frequency, uib, 
drift frequency at 9 = 0, wn, and precession frequency, ui p = (ton)- The trapped 
particle frequencies depend on vn, v±_ and k y . Figure 7(b) shows where in velocity 
space the trapped particle frequencies match the mode frequency for the k y pi=0.6 
mode. The contours indicate that ujd, ujb and lu p ~ O (w) for a thermal particle, and 
that all three resonances lie within the range 0.5vth,e — 3.5vth,e, an d may therefore 
have significant impact. 

This poses several thoughts for analytic theory. Firstly the perturbation changes 
significantly in one bounce period, due to the proximity of u and u>b- Bounce 

% In Ref. [15] a residual instability remained in the absence of magnetic drifts (provided <f> was 
retained), which may be due to a weaker stabilising effect from trapped particles at lower e. 
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averaging, which is often used to simplify the trapped particle response, is therefore 
not appropriate here. Secondly the magnetic drift frequencies are of the same order as 
the mode frequency, and cannot be treated as small. Finally, with all three frequencies 
of the same order as the mode frequency, two or three wave resonances may be 
important. 

4. Conclusions 

Gyrokinetic simulations have found that microtearing modes (MTMs) are unstable in 
both STs and large aspect ratio devices. Recent simulations find that MTMs are also 
unstable in the shallow gradient region just inboard of the MAST H-mode pedestal, 
which may impact on its evolution between ELMs. Analytic theory has proposed 
two different linear drive mechanisms for MTMs: one based on a simple slab model, 
and the other requiring trapped particles. Both mechanisms require a finite rate of 
electron-ion collisions (y e i > 0) for instability. 

A detailed study of the basic linear properties of edge MTMs in MAST has been 
performed using a simplified circular s-a model equilibrium. Consistent with existing 
theories it is found that the mode frequency uj ~ w* e , and the modes are unstable 
only if finite stability thresholds are exceeded in dT e /dr and j3 e . The growth rate 
dependence on u e i, however, is in conflict with existing models. 7 MTM for the dominant 
mode is maximised in the absence of collisions (i.e. at v e i = 0), where the existing 
drive mechanisms should vanish. Trapped particles are essential to drive these MTMs, 
and sensitivity of 7 MTM to the magnetic drift frequency suggests that a drift resonance 
may be involved. 7 MT m rises with the trapped particle fraction, and the k y pi associated 
with the dominant mode drops with increasing s. The mode frequency and the thermal 
trapped electron bounce, precession and drift frequencies are all of the same order. 
To our knowledge this regime has not been addressed by an existing analytic theory. 
In present models the magnetic drifts are typically neglected or treated as small, and 
any trapped particle response is usually obtained using bounce averaging. Neither of 
these approximations are valid here. 

The drive for similar MTMs, at k y pi ~ 0(1), should be enhanced in the high 
magnetic shear region of the edge plateau in tokamaks, and perhaps especially in STs. 
Similar MTMs have also recently been found unstable towards the edge of conventional 
aspect ratio tokamaks including JET [25, 36] and ASDEX Upgrade [17], suggesting 
that this neglected drive mechanism may have wide ranging significance. 
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